Geometric phase of a central spin coupled to an antiferromagnetic environment 
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Using the spin-wave approximation, we study the geometric phase (GP) of a central spin (signal 
qubit) coupled to an antiferromagnetic (AF) environment under the application of an external global 
magnetic field. The external magnetic field affects the GP of the qubit directly and also indirectly 
through its effect on the AF environment. We find that when the applied magnetic field is increased 
to the critical magnetic field point, the AF environment undergoes a spin-flop transition, a first- 
order phase transition, and at the same time the GP of the qubit changes abruptly to zero. This 
sensitive change of the GP of a signal qubit to the parameter change of a many-body environment 
near its critical point may serve as another efficient tool or witness to study the many-body phase 
transition. The influences of the AF environment temperature and crystal anisotropy field on the 
GP are also investigated. 
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The notion of adiabatic geometric phase (GP) of a 
quantum system first discussed in the pioneer work of 
Berry Q was generalized and extended to nonadiabatic 
and noncyclic Q evolutions. It was also generalized 
to mixed states and nonunitary evolutions [4-6]. Re- 
cently, the close relation between the GP's and quantum 
phase transition (QPT) in many-body systems was sug- 
gested and investigated @- In the QPT witness studies, 
an auxiliary spin or signal qubit is introduced to couple 
with a many-body system (environment) and the notion 
of Loschmidt echo (or equivalently the decoherence fac- 
tor of the qubit) [8] or quantum state fidelity Q is used 
to investigate the quantum criticality. The signature of 
QPT is the dramatic decay of the asymptotic value of 
the Loschmidt echo or the quantum state fidelity at the 
critical point. That is, the closer the environment to the 
QPT, the smaller the asymptotic value of the Loschmidt 
echo or the decoherence factor of the signal qubit. Re- 
cently, the GP of an auxiliary qubit induced by a one- 
dimensional XY spin chain (an external environment) 
was calculated and the result was used to study the criti- 
cality of a special case of the transverse Ising model ■ 
It was found that the induced GP in the qubit changes 
dramatically at the critical point [Io| . 

Antiferromagnets subjected to an external magnetic 
field attracted considerable attention over the years [111 ■ 
One of the interesting phenomena in antiferromagnetic 
(AF) materials under an applied magnetic field is the 
magnetic-field-induced spin-flop transition. When the 
applied magnetic field is increased to the critical field 
point, the AF polarization flips into the direction per- 
pendicular to the field. This is called the spin-flop tran- 
sition, a first-order phase transition. The phenomena 
of the spin-flop transition were observed experimentally 



[T3 . fl3j ] . It is thus particularly interesting to investigate 
how a globally applied external field influences the GP of 
a signal qubit that is coupled to an AF environment es- 
pecially when the magnetic field strength is near the crit- 
ical field of the spin-flop transition. It is known that the 
values of the critical magnetic field of the spin- flop transi- 
tion can be obtained using the spin- wave theory. We thus 
apply the spin- wave approximation to deal with the AF 
environment in a globally applied magnetic field at low 
temperatures and low-energy excitations. The spin-wave 
approach allows us to go beyond the Markovian approx- 
imation and the weak-coupling limit in the usual treat- 
ments for system-coupled-to-environment models. The 
influence of the temperature and the crystal anisotropy 
field on the GP will also be investigated. 

The total Hamiltonian of our model [3, [l5[ can be 
written as H = Hs + Hs b + Hb , where 
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are the Hamiltonians of the signal qubit (central spin), 
the coupling, and the environment, respectively. Her g 
is the gyromagnetic factor and \xb is the Bohr magne- 
ton. For simplicity, a significant interaction between the 
central spin and the environment is assumed to be of the 
Ising type with Jo being the coupling constant. J is the 
exchange interaction and is positive for the AF environ- 
ment. We assume that the spin structure of the environ- 
ment may be divided into two interpenetrating sublat- 
tices a and b with the property that all nearest neighbors 
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of an atom on a lie on b and vice versa. Each sublattice 
contains N atoms with spin S, and S a .j (Sf,j) represents 
the spin operator of the ith (jth) atom on sublattice a 
(b). The indices i and j label the AT atoms, whereas the 
vectors S connect atom i or j with its nearest neighbors. 
B represents a uniform external magnetic field applied in 
the z direction. The anisotropy field Ba is assumed to 
be positive, which approximates the effect of the crystal 
anisotropy energy, with the property of tending for posi- 
tive magnetic moment /is to align the spins on sublattice 
a in the positive z direction and the spins on sublattice 
b in the negative z direction. 

Using the Holstein-Primakoff transformation to map 
the spin operators of the AF environment onto bosonic 
operators, considering the situation that the environment 
is in the low-temperature and low-excitation limits, then 
transforming the resultant Hamiltonians to the momen- 
tum space, and finally using the Bogoliubov transforma- 
tion, we obtain in the spin-wave approximation [l5[ 
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where a k (ctu) and (/3k) are the creation (annihilation) 



operators of the two different magnons with wave vector 
k and frequency (o; k respectively, and % = 1 

in Eq. (JS|). For a cubic crystal system in the small k 
approximation 
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where M is the number of nearest neighbors of an atom 
and I is the side length of cubic primitive cell of the sub- 
lattice. 

We assume that the initial density matrix of the to- 
tal systems is separable (i.e., p(0) = IV'KV'I ® Pb, where 
the density matrix of the environment p B = e~ HB / T /Z, 
where Z is the partition function and the Boltzmann con- 
stant has been set to one. The initial state of the central 
spin is described by \4>) — sin(6*o/2)|e) + cos(9 /2)\g) , 
where |e), \g) denote the excited and ground states of 
the qubit and 9q is the polar angle of the initial state 
in the Bloch sphere representation with \g) at the north 
pole. The reduced density matrix of the qubit in the 
thermodynamics limit (i.e., N — > oo) can be obtained, 
following the calculation in Ref. [l5| , to be 
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The decoherence time tq is given by [If 



V27r/[W(»? + +T)]. where 
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with x = kl. The decoherence time depends on the cou- 
pling strength to the environment, on the structure pa- 
rameters and temperature of the environment, and on 
the external magnetic field. 

To calculate the GP for the qubit undergoing nonuni- 
tary evolution, we use the gauge invariant expression de- 
rived in [j| 
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where r denotes the total evolution time, and £fc(r) and 
\4>k{T~)) are the eigenvalues and corresponding eigenvec- 
tors of the reduced density matrix pit) . One may view 
the GP factor defined in Eq. ([9]) as a weighted sum over 
the phase factors pertaining to the eigenvectors of the 
reduced density matrix. The eigenvalues and eigenstates 
of p{t) in Eq. (J7J can be obtained as 
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Substituting Eqs. jlO]) and (JTTJ into Eq. © and carry- 
ing out some simple manipulations, we obtain the GP 
acquired at the time r = 27r / \gpsB) of a cyclic period 
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The AF environment causes decoherence and influences 
the GP of the signal qubit. Our results obtained using the 
spin- wave approximation are valid only for not too large 
temperatures, much smaller than the Neel temperature. 
According to neutron diffraction studies, the Neel tem- 
perature of antiferromagnet TbAuIn is 35 K (l6l | . Some 
antiferromagnets may have higher Neel temperatures. So 
in the following analysis, the environmental temperature 
is restricted below T/gp B = 2.5Tesla (i.e., T w 3.4K). 
, We, next, discuss the GP acquired at the time r = 
) 2-7T / 'igpsB) after the qubit completes a cyclic evolution 
/ when it is isolated from the AF environment. If the qubit 
(7)s really isolated from the AF environment, i.e., Jo = 0, 
then the decoherence time r n — > oo. In this case, we ob- 
tain from Eq. (|12[) the well-known accumulated geometric 
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FIG. 1: (Color online) The GP of the signal qubit as a 
function of the polar angle 8o for Jo = (solid curve), 
Jo = 2. 5 J, T/g/j,B = 0.8Tesla (dashed curve), and Jo = 2.5J, 
T/gfiB = 1.2Tesla (dot-dashed curve). Other parameters 
are M = 6, MJ/gfj, B = 40Tesla, B A = O.lOTesla, and 
B = 0.5 Tesla 



phase of a spin-1/2 particle processing around an exter- 
nal magnetic field in each cycle as: $ = tt(1 — cos#o). If 
6*o = f, we have $ = n. 

Figure [1] shows the GP of the qubit as a function 
of the initial polar angle 8 Q for different temperatures. 
The external magnetic field is chosen as B = 0.5 Tesla 
and other parameters are with Ba = 0.1 Tesla and 
MJ/gfiB = 40 Tesla. The solid curve corresponds to the 
qubit isolated from the AF environment (Jo = 0) [i.e., 
$ = 7r(l — cos 6>o)] . Comparing with the isolated case, we 
can see from Fig. [1] that the influence of the environment 
makes the GP curve deviate steeply from the isolated 
case in the region around do = tt/2. The deviation of the 
GP is positive for 6q > tt/2 and is negative for 9o < tt/2. 
Furthermore, the higher the temperature is, the steeper 
the curve is. For three special initial states with polar 
angles do = 0, 0q = § j an d = the GP's are not 
affected by the AF environment. 

In Fig. [21 we plot the GP as a function of external mag- 
netic field B for different temperatures with 6q = 1.3. 
The solid curve corresponds to the isolated case ( J = 0) 
[i.e., $ = 7r(l — cos 1.3) = 0.737r], which is independent 
of B for B ^ 0. It is obvious from Fig. [5] that the curve 
with higher temperature shows greater deviation from 
$ = 0.737T. For a given strength of the external mag- 
netic field, the decoherence time To of the off-diagonal el- 
ements of the reduced density matrix decreases with the 
increase of temperature. This reduces the GP obtained 
by the qubit @. The influence of the external magnetic 
field on the GP (TT21 can be categorized into two com- 
peting ways: one through the decoherence time tq and 



the other one through the quasiperiod r = 2tt/ (gfisB). 
If the external magnetic field increases, the decoherence 
time To decreases Then it suggests that the influ- 

ence of the environment on the GP increases with the 
increasing magnetic field. On the other hand, with the 
increase of the external magnetic field B, the quasiperiod 
t decreases. It then takes less time for the environment 
to influence the GP in a quasicycle. In this sense, the 
influence of environment on the GP decreases with the 
increase of the external magnetic field. Hence, when the 
external magnetic field is very small, the decoherence 
time t is finite [13, EH but the qubit needs a rather 
long time t > To to complete a quasicycle. Since the 
integrand with the factor exp(— 2t 2 /t^) in Eq. (fl~2|) be- 
comes very small for time t considerably larger than To, 
the integration over t from To to r thus contributes very 
little to the whole GP of the signal qubit as compared 
to the isolated case. As a result, the dashed and dot- 
dashed curves at very small B show large deviation from 
the value $ = 0.737T of the isolated case. One can fur- 
thermore obtain from Eq. (fT2|) that the GP near B = 



behaves linearly in B as lim $ = c\B 1 where c\ is in- 
dependent of B. As one increases the strength of the 
magnetic filed, the decoherence time decreases but the 
quasiperiod also decrease. Usually, the second effect is 
more important in typical values of the external mag- 
netic field. Therefore, the low-temperature dashed curve 
approaches $ = 0.737T with the increase of the external 
magnetic field. When the strength of the magnetic field 
is increased further toward the critical point where the 
spin-flop transition occurs, B = B c m 2.83 Tesla for the 
parameters used in Figs. [5] and [31 the decoherence time 
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FIG. 2: (Color online) The GP of the signal qubit as a func- 
tion of the external magnetic field B for Jo = (solid curve), 
Jo = 2. 5 J, T/g/j,B = 0.8 Tesla (dashed curve), and Jo = 2. 5 J, 
T/gfiB ~ 1.5 Tesla (dot-dashed curve). Other parameters are 
0o = 1.3, M = 6, MJ/g^B = 40 Tesla, and B A = O.lOTesla. 
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approaches zero abruptly [15| . An analytical expression 
for the decoherence time near the critical field can be 
found in Eq. (39) of Ref. [l5[. We note that this result us- 
ing spin- wave theory is valid for B < B c (i.e., before the 
spin- flop transition). Thus when B — > B c from below, the 
environment exerts a great influence on the qubit so that 
its coherence is destroyed completely, but the quasiperiod 
is still finite. Such a great influence thus also appears on 
the GP and causes a sudden deviation from $ = 0.737T for 
the dashed curves in Fig. [5] when the external magnetic 
field approaches B c — 2.83 Tesla. The drastic change in 
the induced GP of the qubit at the critical point is due 
to the sensitivity of the many-body (AF) environment to 
the parameter (the magnetic field) change near its crit- 
ical point. Furthermore, from Eq. (39) of Ref. [l5[ and 
Eq. (H2), one can obtain the GP near the critical field 
scaling as lim $ = C2(B C — B) 1 ^ 4 . Here C2 is in- 

B^B c (B<B a ) 

dependent of the external magnetic field B. The linear 
behaviors of the GP near B = and the abrupt changes 
near B = B c can be clearly observed in Fig. [3J 

In Fig. |31 the GP of the qubit (central spin) is plotted 
as a function of B and T for the value of anisotropic field 
Ba = 0.10 Tesla. Obviously, the curves in Fig. [5] arc 
cross-section curves of the curved surface in Fig. at dif- 
ferent temperatures. As mentioned previously, one can 
observe from Fig. [3] that the GP decreases from the value 
<J> = 0.737T with the increase of the temperature. Besides, 
for X '/ 'gfiB > 0.5 Tesla, it decreases quickly when B ap- 
proaches zero. As the external magnetic field approaches 
the critical point of B c ss 2.83 Tesla, the GP changes 
drastically to zero. Figure S] is similar to Fig. [31 except 
that the anisotropy field is at Ba = 0.15 Tesla. Then 
the spin-flop transition occurs at B = B c sa 3.47 Tesla. 
In such a case, the GP does not show abrupt change 
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FIG. 3: (Color online) The GP of the signal qubit versus the 
external magnetic field B and the environment temperature 
T. Other parameters are 9o = 1.3, M = 6, MJ/gfiB = 
40 Tesla, J = 2. 5 J, and B A = 0.10 Tesla. 
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FIG. 4: (Color online) Same as Fig. 3, except Ba = 
0.15 Tesla. 



even when B approaches 3 Tesla. Comparing Fig. [3] with 
Fig. [4] we find that the large crystal anisotropy field sup- 
presses the influence of the AF environment on the GP. 

In conclusion, we investigated the influence of an AF 
environment on the GP of a signal qubit (central spin) . 
Such an influence is enhanced by increasing the temper- 
ature or decreasing the crystal anisotropy field. The de- 
pendence of the GP on the external magnetic field B 
involves two competing factors: quasiperiod r and deco- 
herence time To- A larger B implies a smaller t which 
indicates that the GP acquired at time r by the qubit is 
less influenced by the environment in terms of the dura- 
tion of interaction time with the environment. But the 
larger B also implies smaller tq which indicates that the 
decays in the integrand of the GP (|12p become greater. 
When B approaches the region near the critical field of 
the spin-flop transition, To — > but r is finite. The qubit 
acquires no GP with To — > and thus the GP changes 
abruptly at the critical field of the spin-flop transition. 
This sensitive change of the GP of a signal qubit to the 
parameter change of a many-body environment near its 
critical point may serve as, in addition to the Loschmidt 
echo and quantum state fidelity, another efficient tool or 
witness to study the phase transition. 
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